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ABSTRACT 


Wring  the  shuttle's  Spacelab-2  mission,  the  Un4_uer#-ity  of  Iowa-'* 
Plasma  Diagnostics  Package  (PDP)  was  released  from  the  shuttle  to 
freefly.  At  times  during  this  freeflight  when  the  PDP  was  magnetic- 

V- 

ally  connected  to  the  shuttle,  -Stanford'^  Fast  Pulsed  Electron  Genera¬ 
tor  (FPEG),  located  in  the  shuttle  cargo  bay,  ejected  a  1  keV  -  50  mA 
electron  beam.  The  plasma  wave  instrument  on  board  the  PDP  detected 
intense  whistler-mode  radiation  during  these  beam  ejections.  This 
paper  presents  a  study  of  a  whistler-mode  emission  detected  during  one 
particular  continuous  electron  beam  firing.  Calculations  indicate 
that  the  beam  radiated  approximately  1.6  mW  in  the  whistler  mode  as 
the  beam  traversed  the  200  meters  from  the  shuttle  to  the  PDP.  The 
emissivity  also  decreased  by  about  a  factor  of  10  over  this  same  dis- 

7  '  '  7"'  / 

tance.  The  measured  wave  powers  are  10'  greater  than  wave  powers 
expected  from  incoherent  Cerenkov  radiation,  verifying  that  the  radia¬ 
tion  is  generated  by  a  coherent  process.  Estimates  of  the  emissivity 
based  on  measured  electric  field  intensities  in  the  beam  indicate  that 
the  whistler-mode  noise  is  produced  by  radiation  from  electron  bunches 
created  by  an  electrostatic  beam-plasma  instability. 


.  i . 

N 


INTRODUCTION 


In  this  paper  we  describe  an  electron  beam  experiment  performed 
on  the  Spacelab-2  mission.  The  Spacelab-2  flight,  which  was  launched 
on  July  29,  1985,  included  an  electron  accelerator  called  the  Fast 
Pulsed  Electron  Generator  (FPEG)  from  Stanford  University,  and  a 
spacecraft  called  the  Plasma  Diagnostics  Package  (PDP)  from  the 
University  of  Iowa.  During  a  6-hour  period  on  August  1,  1986,  the  PDP 
was  released  from  the  shuttle  to  investigate  plasma  effects  in  the 
vicinity  of  the  shuttle.  During  the  PDP  free  flight,  the  shuttle  was 
maneuvered  so  that  the  PDP  passed  near  magnetic  field  lines  connected 
to  the  shuttle.  Four  such  magnetic  conjunctions  were  achieved.  Dur¬ 
ing  one  of  these  magnetic  conjunctions  a  1  keV  -  50  mA  electron  beam 
was  continuously  ejected  from  the  shuttle  so  that  radiation  effects 
could  be  monitored  as  the  PDP  passed  near  the  magnetic  field  line 
carrying  the  beam.  Figure  1  shows  a  frequency  vs.  time  spectrogram 
from  the  PDP  plasma  wave  instrument  during  this  electron  beam  event. 
The  funnel-shaped  signal  extending  from  the  electron  cyclotron  fre¬ 
quency,  fc,  down  to  approximately  30  kHz  is  whistler-mode  radiation 
from  the  electron  beam.  This  whistler-mode  radiation  was  first 
described  by  Gurnett  et  al.  [1986].  Our  objectives  in  this  paper  will 
be  to  determine  the  total  radiated  power  from  the  beam  and  compare  it 
with  the  power  predicted  by  various  whistler-mode  radiation  mecha¬ 
nisms. 
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Whlstler-mode  emissions  are  known  to  be  produced  by  both  artifi¬ 
cial  and  natural  electron  beams.  Some  of  the  early  observations  of 
emissions  from  artificial  electron  beams  were  made  in  the  early  seven¬ 
ties  by  the  University  of  Minnesota  group  using  the  Echo  rocket  experi¬ 
ments  [Cartwright  and  Kellogg,  1974;  Kellogg  et  al.,  1976;  Monson  et 
al.,  1976;  Winckler,  1980].  Observations  of  beam-generated  emissions 
were  also  made  on  the  joint  Franco-Soviet  ARAKS  beam  experiments 
[Dechambre  et  al.,  1980].  The  first  space  shuttle  electron  beam  exper¬ 
iment  was  performed  in  March,  1982,  as  part  of  the  STS-3  mission.  On 
this  flight,  the  PDP  was  carried  on  the  shuttle  remote  manipulator  arm 
and  the  FPEG  was  used  to  produce  an  artificial  electron  beam.  During 
electron  gun  firings,  a  possible  whistler-mode  signal  was  detected  from 
the  electron  beam  [Shawhan  et  al. ,  1984].  In  December  of  1983,  the 
Phenomenon  Induced  by  Qiarged  Particle  JJeams  (PICPAB)  experiment  was 
carried  as  part  of  the  Spacelab-1  mission.  The  PICPAB  experiment  was 
specifically  designed  to  investigate  radiation  from  electron  beams 
injected  in  the  ionosphere.  Like  the  other  previous  experiments, 
whistler-mode  signals  were  again  detected  when  the  electron  gun  was 
fired  [Beghin  et  al.,  1984], 

Whistler-mode  radiation  is  also  produced  in  the  auroral  zone  in 
association  with  the  f ield-al igned  electron  beams  that  are  responsible 
for  the  aurora  [Gurnett,  1966].  This  radiation  is  usually  called 
auroral  hiss.  Both  upward  and  downward  propagating  auroral  hiss  has 
been  observed  [Mosier  and  Gurnett,  1969].  The  downward  propagating 
auroral  hiss  is  associated  with  downward  moving  electron  beams  with 
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characteristic  energies  of  a  few  hundred  eV  [Gurnett,  1966;  Hartz, 

1969;  Gurnett  and  Frank,  1972;  Laaspere  and  Hoffman,  1976].  The  upward 
propagating  auroral  hiss  often  has  a  V-shaped  spectrum  called  a 
"saucer"  (Smith,  1969;  Mosier  and  Gurnett,  1969;  James,  1976]  or  a 
"funnel"  [Gurnett  et  al.,  1983].  Upward  propagating  auroral  hiss  has 
been  observed  in  association  with  upward  moving  field-aligned  electron 
beams  [Lin  et  al.,  1984].  The  characteristic  frequency-time  shape  of 
the  "saucer"  or  "funnel"  is  a  propagation  effect  that  occurs  for 
whistler-mode  waves  propagating  near  the  resonance  cone.  Gurnett  et 
al.  [1986],  in  their  Initial  analysis  of  the  Spacelab-2  results,  showed 
that  a  similar  wave  propagation  effect  was  observed  for  the  whistler¬ 
mode  signals  generated  by  the  SL-2  electron  beam. 


i 
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ELECTRIC  FIELD  POLARIZATION 


To  confirm  that  the  radiation  from  the  Spacelab  2  electron  beam  is 
propagating  in  the  whistler  mode,  we  compare  the  electric  field  polar¬ 
ization  with  the  polarization  expected  for  the  whistler  mode.  The 
whistler  mode  has  a  polarization  that  depends  on  the  wave  frequency,  f, 
the  wave  normal  angle,  0,  the  cyclotron  frequency,  fc,  and  the  plasma 
frequency,  fp.  Using  cold  plasma  theory  [Stix,  1962],  the  electric- 
field  and  index  of  refraction  vectors  can  be  calculated  as  a  function 
of  f,  0,  fc,  and  fp.  For  the  whistler  mode  the  index  of  refraction 
varies  strongly  as  a  function  of  the  wave  normal  direction.  This  vari¬ 
ation  can  be  represented  by  an  index  of  refraction  surface,  n(0),  which 
is  a  surface  that  is  defined  by  the  loci  of  index  of  refraction  vectors 
at  different  wave  normal  angles.  Figure  2  shows  the  index  of  refrac¬ 
tion  surface  for  the  whistler  mode.  At  a  limiting  wave  normal  angle, 
known  as  the  resonance  cone  angle,  0Res ,  the  index  of  refraction  goes 
to  infinity.  The  resonance  cone  angle  is  given  by  tanz0Res  =  -P/S, 
where  P  =  1  -  fp2/f2  and  S  =  1  -  fp2/(f2  -  fc2).  As  the  wave  normal 
approaches  the  resonance  cone,  the  electric  field  E  becomes  linearly 
polarized  with  E  parallel  to  n.  In  this  limit  the  electric  field  is 
quasi-eiectrostatic  and  the  group  velocity,  vg,  is  perpendicular  to  ~ 
n  (see  Figure  2). 

In  a  previous  paper  [Curnett  et  al.,  1986],  the  funnel-shaped 
frequency  versus  time  pattern  of  the  radiation  from  the  SL-2  electron 
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beam  was  explained  as  a  frequency  dependent  propagation  effect  for 
whistler-mode  emissions  from  the  electron  beam.  Figure  3  is  a  sketch 
demonstrating  how  this  effect  works.  Consider,  first,  a  point  source 
of  whistler-mode  radiation,  with  the  radiation  propagating  near  the 
resonance  cone.  As  the  wave  frequency  increases,  the  resonance  cone 
angle  decreases  and  the  ray  path  direction,  Vg,  becomes  increasingly 
oblique  to  the  magnetic  field,  approaching  90°  as  the  frequency  approa¬ 
ches  the  electron  cyclotron  frequency.  Typical  ray  paths  for  a 
whistler-mode  emission  propagating  near  the  resonance  cone  are  shown  in 
the  figure.  As  a  spacecraft  approaches  the  source,  emissions  near  the 
gyrof requency  are  detected  first,  since  their  ray  paths  are  almost  per¬ 
pendicular  to  the  beam.  Lower  and  lower  frequencies  are  then  detected 
as  the  distance  between  the  spacecraft  and  source  decreases.  This 
frequency  dependent  propagation  effect  creates  a  funnel-shaped  emission 
pattern  on  a  frequency  versus  time  spectrogram.  By  comparing  a  modeled 
emission  pattern  using  cold  plasma  theory  to  the  funnel-shaped  pattern 
shown  in  Figure  1,  Gurnett  et  al .  [1986]  have  confirmed  that  the 
whistler-mode  emission  from  the  SL-2  electron  beam  is  propagating  up¬ 
ward  from  the  shuttle  near  the  resonance  cone.  The  distinct  funnel 
boundary  is  defined  by  the  waves  that  propagate  directly  from  the  beam 
origin,  the  shuttle.  The  funnel  on  the  spectrogram  appears  filled-in 
since  the  radiation  is  being  emitted  by  a  line  source  extending  upward 
from  the  shuttle  rather  than  a  point  source.  Since  the  waves  are 
upward  propagating,  k*vjj  >  0,  where  vjj  is  the  beam  velocity  vector. 

In  order  to  provide  further  confirmation  that  the  radiation  from 
the  SL-2  electron  beam  is  propagating  near  the  resonance  cone,  an 


additional  test  was  performed.  This  test  compares  model  electric- 
field  directions  in  the  PDP  spin  plane  to  their  actual  directions  as 
measured  by  the  PDP  plasma  wave  instrument.  To  perform  this  test  a 
computer  program  was  developed  that  calculates  the  angle,  <J> ,  between 
the  projection  of  a  model  electric  field  onto  the  spin  plane  and  a 
fixed  reference  direction.  The  fix  reference  direction  selected  was 
the  spin  plane  projection  of  the  spacecraft-sun  vector.  To  compute  b, 
the  group  velocity  was  assumed  to  be  directed  from  a  point  on  the 
beam  toward  the  PDP  with  the  electric  field  vector,  E,  at  an  angle 
0Res  relative  to  the  beam.  This  field  geometry  is  the  expected  con¬ 
figuration  for  an  upward  propagating  whistler-mode  wave  near  the 
resonance  cone.  Figure  4  shows  the  corresponding  geometry  of  E,  vg 
and  Tc. 

The  electric-field  directions  in  the  spin  plane  calculated  using 
the  model  described  above  are  compared  to  the  measured  electric-field 
directions  found  from  spin  modulation  maximums  in  the  receiver  data. 
The  spin  modulation  maximums  occur  when  the  PDP  electric  antennas  are 
aligned  with  the  measured  electric  field  in  the  spin  plane,  thus 
allowing  a  direct  determination  of  this  measured  electric  field  direc¬ 
tion.  Figure  5  shows  the  results  of  this  comparison  at  four  frequen¬ 
cies:  562,  311,  178,  and  100  kHz.  This  figure  shows  the  phase  angle, 

t»  between  the  projected  electric  field  and  the  sun  vector  as  a  func¬ 
tion  of  time.  The  dots  represent  the  modeled  electric-field  direc¬ 
tions  computed  assuming  a  resonance  cone  propagation  scheme  while  the 
X's  represent  the  measured  electric-field  directions.  The  close 
agreement  between  the  computed  and  measured  electric  field  directions 
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provides  strong  confirmation  that  the  waves  are  propagating  near  the 
resonance  cone  and  in  the  beam  direction  (i.e.,  k'V^  >0),  as 
indicated  in  Figure  4. 
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EMITTED  POWER 


In  this  section  we  estimate  the  total  power  radiated  from  the 
beam  in  the  whistler  mode.  By  comparing  the  radiated  power  to  the 
total  power  in  the  beam,  the  efficiency  of  the  wave-beam  interaction 
can  be  estimated  and  compared  with  various  generation  mechanisms. 

The  power  emitted  from  the  beam  in  the  whistler  mode  can  be  estimated 
by  integrating  the  Poynting  flux  over  a  surface  surrounding  the  beam. 
An  inherent  difficulty  with  this  calculation  is  the  determination  of 
the  phase  and  magnitude  of  the  electric  and  magnetic  fields  in  the 
Poynting  flux,  S  =  ExH.  Since  three  axis  measurements  are  not  avail¬ 
able  and  since  phase  measurements  were  not  made  the  Poynting  vector 
cannot  be  determined  directly.  Computing  the  Poynting  flux  from  the 
electric  field  alone  is  complicated  by  the  fact  that  near  the  reso¬ 
nance  cone  the  whistler  mode  is  quas i-e iectrostat ic  and  the  ratio  of 
the  electromagnetic  to  electrostatic  components  is  an  exceedingly 
sensitive  function  of  the  wave  normal  angle.  To  compute  the  wave 
normal  angle  we  assume  that  the  radiation  is  produced  by  the  Landau 
resonance,  i.e.  ,  <n /k ;f  =  v^.  Since  the  beam  velocity  is  known,  this 
condition  gives  a  well-defined  value  for  the  wave  normal  direction. 
The  fact  that  the  radiation  is  propagating  in  the  same  direction  as 
the  beam  (k*V}j  >  0)  provides  a  strong  indication  that  the  Landau 
resonance  is  involved.  For  example,  the  s  =  -1  cyclotron  resonance 
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produces  radiation  propagating  in  the  opposite  direction  of  the  beam. 
Also,  as  will  be  discussed  later,  the  Landau  resonance  gives  the  best 
agreement  with  the  measured  electric  to  magnetic  field  ratios. 

To  carry  out  the  Poynting  vector  calculation  we  note  that  E  lies 
almost  entirely  in  the  plane  defined  by  n  and  the  geomagnetic  field 
(see  Figure  4).  The  electrostatic  and  electromagnetic  components  of  E 
are  then  given  by  Ec  cos  A0  and  E0  sin  A8 ,  respectively,  where  A9  is 
the  angle  between  E  and  n,  and  E0  is  amplitude  of  the  total  electric 
field.  The  angle  A0  is  determined  by  the  Landau  resonance  condition 
and  cold  plasma  theory.  The  Landau  resonance  condition  specifies  the 
component  of  n  parallel  to  the  geomagentic  field,  i.e., 


nn  =  n  cos  0  =  c/vb 


where  c  is  the  speed  of  light.  For  a  1  keV  electron  beam  moving 
parallel  to  the  magnetic  field  n u  is  approximately  15. 9.  A  program  was 
written  that  solves  Equation  (1-20)  of  Stix  [1962]  for  the  magnitude 
and  directions  of  n  and  F.  Equation  (1-20)  of  Stix  is 


S  -  n2cos20 
i  D 

n2cos  0  sin  0 


-i  0  n^cos  9  sin  0 

S  -  n2  0 

0  P  -  n2sin20 


=  0  (3) 
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P  =  1  -  E  -EJ~ 
k  u 


(8) 


The  quantities  u>pk  and  (ocr  are  the  plasma  and  cyclotron  frequencies 
for  species  k  while  ek  is  the  sign  of  the  charge  of  species  k.  Using 
this  program,  n  and  A0  at  a  particular  wave  frequency  can  be  calculated 
consistent  with  nn .  Since  A0  Is  now  determined,  the  electrostatic  and 
electromagnetic  components  of  F,  can  be  calculated.  The  calculated  A0 
values  are  small,  typically  ranging  from  .06°  to  1.1°  from  31.1  kHz  to 
562  kHz,  and  indicates  that  the  wave  is  nearly  electrostatic.  It  is 
easy  to  show  that  the  magnitude  of  the  Poynting  vector  is  given  by 


n 


1/2 


S  = 


(_£)  (A2  +  B2)l/2 

2  un 


(9) 


where  A  =  1  -  cos2A0  and  R  =  sin  A0  cos  AO.  In  the  derivation  of 
Equation  (9)  Faraday's  Law,  n  x  F,  =  cB,  was  used  to  eliminate  the 
magnetic  field  in  the  E  x  H  term.  Consequently,  |s|  is  only  dependent 
on  n  and  E  which  are  we L 1  defined  quantities.  Note,  that  as  0 


approaches  the  resonance  cone  angle,  n  and  F  become  parallel  and  |s| 
goes  to  zero.  This  behavior  near  0Res  is  similar  to  an  expression 
derived  by  Mosier  and  Gurnett  [1*571]  in  their  paper  addressing  Poynting 
flux  measurements  of  VLF  hiss  emissions. 

Figure  6  shows  the  PDP  trajectory  relative  to  the  beam.  As  can 
be  seen,  the  PDP  trajectory  was  nearly  perpendicular  to  the  beam.  At 
closest  approach,  it  passed  within  about  3  meters  of  the  beam  at  a 
distance  of  roughly  200  meters  from  the  shuttle.  To  compute  the  total 
radiated  power,  the  Poynting  flux  is  integrated  over  an  imaginary  sur¬ 
face  perpendicular  to  the  beam  that  includes  the  PDP  trajectory. 
Assuming  that  the  sampled  intensities  along  this  trajectory  are  ''on- 
stant  around  an  annular  ring  of  the  area,  dA,  centered  on  the  beam,  the 
radiated  power  from  the  beam  segment  can  be  obtained  by  evaluating  the 
integral  p  =  fS||2irRdR,  where  S  n  is  the  field-aligned  component  of  the 
Poynting  vector  and  R  is  the  perpendicular  distance  from  the  beam  to 
the  PDP.  Note  that  the  evaluation  of  this  integral  will  vield  two 
values  for  the  radiated  power:  one  value  from  the  inbound  pass  where 
the  limits  of  integration  extend  from  R  =  «  to  R  =  0  and  one  value  from 
the  outbound  pass  where  the  limits  of  the  integration  now  extend  from  R 
=  0  to  R  =  -°°.  Figure  7  shows  the  average  power  spectral  density  from 
these  two  passes  as  a  function  of  frequency.  The  error  bars  shown  in 
the  figure  represent  the  standard  deviations  of  the  power  values.  Note 
that  the  power  spectral  density,  dP/df,  is  on  the  order  of  10-9  W/Hz  in 
the  frequency  range  extending  from  30  kHz  to  1  MHz.  The  total  emitted 
power  in  the  200  meter  beam  segment  from  the  shuttle  to  the  PDP  is 
found  to  be  P  =  1.6  mW.  Assuming  that  the  power  is  emitted  uniformly 


along  the  beam,  the  radiated  power  per  unit  length,  dP/dH,  is  approx¬ 
imately  8  x  10-6  W/m.  Since  the  total  power  of  the  beam  was  50  W,  the 
beam  converted  approximately  3.2xlO-5  of  its  power  to  whistler-mode 
radiation  in  the  first  200  meters.  If  the  beam  dissipated  energy  only 
via  the  whistler-mode,  the  beam  would  only  propagate  approximately  6000 
km  before  all  of  the  beam  energy  would  be  converted  to  radiation. 

The  linear  emissivitv  of  the  whistler-mode  radiation,  dP/dfdH, 
from  different  locations  along  the  beam  can  also  be  calculated.  To 
calculate  the  linear  emissivity,  a  knowledge  of  a  signal's  exact  source 
location  from  the  beam  is  required;  however,  by  using  the  ray  path,  the 
source  of  the  signal  at  a  particular  point  along  the  PDP  trajectory  can 
be  located.  The  power  radiated  from  an  infinitesimal  beam  radiation 
source,  Al ,  is  P  =  /  2irRdH,  where  Sj_  is  the  perpendicular  component 
of  the  Poynting  vector  measured  at  the  perpendicular  distance  R  from 
the  beam  and  corresponds  to  the  Poynting  flux  emitted  from  a  cylinder 
of  radius,  R,  and  length,  d£,  surrounding  the  beam.  The  linear  emis¬ 
sivity  from  this  source,  dP/dfdi,  is  then  obtained  by  using  the  dif¬ 
ferential  form  of  the  power  intergal.  The  calculated  linear  emissivitv 
of  the  whistler-mode  waves  is  shown  in  Figure  8.  Note  that  the  emis¬ 
sivitv  drops  by  a  factor  of  ten  from  100  ro  200  meters  along  the  beam. 
This  decrease  in  emissivitv  indicates  that  the  efficiencv  of  whistler¬ 
mode  generation  decreases  with  increasing  distance  along  the  beam  and 
that  the  generation  mechanism  is  capable  of  dvnami c  changes  in  tens  of 
meters.  If  the  emissivitv  continues  to  drop  at  the  rate  observed 
between  100  to  200  meters,  the  radiation  would  be  undetectable  v^v  the 


PDP  at  source  distances  more  than  about  1  km  from  tN>  shuttle. 


This 


result  may  explain  why  DE-1,  which  was  magnetically  connected  to  the 


shuttle  during  a  gun  firing  on  the  STS-3  mission,  did  not  see  beam¬ 
generated  whistler-mode  radiation  in  the  vicinity  of  the  streaming 
electrons  [ Inan  et  al.,  1984],  From  the  SL-2  measurements,  it  appears 
that  strong  whistler-mode  emissions  are  probably  generated  only  in 
close  proximity  to  the  source  of  the  beam. 

As  mentioned  earlier,  the  electric  and  magnetic  field  measurements 
also  provide  further  evidence  that  the  whistler-mode  waves  are  gen¬ 
erated  via  a  Landau  resonance  process.  This  argument  involves  a  com¬ 
parison  of  computed  and  measured  cB/E  ratios.  Assuming  a  specific  res¬ 
onance  condition  and  using  the  solution  of  Equations  (1-20)  of  Stix 
[1962],  a  unique  value  for  n  and  A0  can  be  computed.  Faraday's  law  can 
then  be  used  to  obtain  the  relationship 

n  x  E  =  cB  (10) 

where  E  is  in  the  electric  component  and  B  is  in  the  magnetic  compo¬ 
nent  of  the  whistler-mode  waves.  For  the  assumed  field  geometry, 
Equation  (10)  can  be  rewritten  as 

n  E0  Sin  A0  =  cBc  , 


or 


n  Sin  A0 


(ID 


Ising  Equation  (11),  we  can  compute  n  sin  A0  for  various  resonance  con¬ 
ditions  and  compare  this  ratio  with  the  measured  cB/E  ratio.  The 
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spectrum  analyzer  used  with  the  PDP  search  coil  can  only  provide  mea¬ 
surements  up  to  178  kHz;  therefore,  the  magnetic  to  electric  field 
ratio  can  only  be  obtained  in  the  66  kHz,  100  kHz,  and  178  kHz  frequency 
channels.  Also,  the  measured  values  of  B  at  high  frequencies  using  the 
search  coll  are  highly  uncertain,  due  to  inaccuracies  in  the  calibration 
of  the  Instrument.  The  preflight  calibration  was  performed  by  placing  a 
calibration  coll  in  the  search  coil  and  surrounding  the  system  in  a  u~ 
metal  can.  A  problem  arises  at  high  frequencies  (>10  kHz),  where  fre¬ 
quency  dependent  capacitences  and  inductances  affect  the  current  and  the 
expected  value  of  B  from  the  calibration  coils.  Unfortunately,  post¬ 
flight  calibrations  under  more  ideal  condition  (specifically,  without 
the  p-metal  can)  have  failed  to  reproduce  the  preflight  calibrations. 
This  suggests  that  the  high  frequency  gain  of  the  search  coil  may  have 
shifted  during  the  flight.  Our  current  best  estimates  are  that  B  (and 
cB/E)  are  accurate  only  to  within  a  factor  of  2  -  4  at  high  frequencies. 
The  range  of  measured  cB/E  values  lies  between  1.1  and  15.1.  Assuming  a 
Landau  resonance,  n  sin  A0  is  computed  to  be  .54,  .52,  and  .54  for  56 
kHz,  100  kHz,  and  178  kHz,  respectively.  Note  that  these  values  lie 
just  outside  the  range  of  measured  cB/E  values,  and  fall  in  the  range 
when  considering  the  factor  of  2  -  4  uncertainty  in  the  calibrations. 

For  an  s  =  +1  cyclotron  resonance,  however,  n  sin  A0  is  computed  to  he 
between  .05  to  .OP,  for  56  kHz,  1  Of)  kHz,  and  178  kHz.  These  values  are 
about  a  factor  of  20  smaller  than  the  lowest  measured  oB/E  value. 

Similar  computed  values  are  obtained  lor  the  s  =  -l  cvclotron  resonance. 
These'  comparisons  show  that  the  measured  cB/E  rat  io  is 
expected  for  a  I.nndau  resonance. 


closest  to  those 


INCOHERENT  GENERATION  MECHANISMS 


From  the  power  measurements  alone  it  is  not  clear  whether  the  beam¬ 
generated  whistler-mode  radiation  detected  by  the  PDP  during  the  SL-2 
mission  results  from  a  coherent  or  incoherent  generation  process.  A 
coherent  mechanism  involves  large  numbers  of  particles  acting  together 
to  generate  the  emitted  waves.  The  total  power  from  a  coherent  source 
goes  as  N2  ,  where  N  is  the  number  of  particles  in  coherence.  Common 
coherent  sources  are  plasma  instabilities,  lasers  and  radio  antennas. 
Incoherent  mechanisms  involve  particles  that  are  radiating  independent¬ 
ly.  The  power  from  the  individual  radiators  must  be  added  to  get  the 
total  power  emitted;  thus  the  total  power  is  proportional  to  N,  the 
number  of  radiators.  A  common  incoherent  source  is  an  incandescent 
light  bulb.  This  section  and  the  next  describe  possible  incoherent  and 
coherent  mechanisms  for  generating  whistler-mode  radiation. 

One  possible  incoherent  mechanism  involves  Cerenkov  radiation. 
Cerenkov  radiation  is  caused  by  charged  particles  moving  with  velocities 
greater  than  the  phase  speed  of  the  wave  in  the  medium.  The  whistler¬ 
mode  waves  from  the  SL-2  electron  beam  are  propagating  near  the  reso¬ 
nance  cone  with  large  index  of  refractions,  typically  n  ~  30  to  500. 

The  phase  speed  of  the  wave  is  then  much  less  than  the  speed  of  the  1 
keV  electron  beam.  Since  the  beam  electrons  are  moving  faster  than  the 
phase  velocity  of  the  whistler  mode,  Cerenkov  radiation  should  be 


18 

The  measured  whistler-mode  power  from  the  beam  will  now  be  compared 

to  the  calculated  power  from  Cerenkov-radiation ,  assuming  that  the  beam 

electrons  are  incoherent  radiators.  Our  calculations  are  similar  to 

those  performed  by  Jorgenson  [1968]  and  Taylor  and  Shawhan  [1973],  who 

both  calculated  the  power  from  this  process  and  compared  it  to  the 

radiated  powers  from  VLF  hiss.  Mansfield  [1967]  derived  an  equation 

that  gives  the  power  spectral  density  radiated  from  a  single  electron 

moving  through  an  ambient  ionized  gas  with  a  speed  greater  than  the  wave 

phase  speed.  For  an  incoherent  mechanism,  the  total  power  radiated  from 

dP 

the  beam  is  the  power  radiated  from  each  electron  (-n-r)  ,  added  up  over 

e 

all  the  electrons  in  a  given  beam  volume,  Nv:  (4^)  =  n(4?-)  • 

dt  total  of  e 

Using  Mansfield's  formula,  the  radiated  power  from  each  beam  electron 

can  be  calculated  and  Is  shown  in  Figure  9.  In  obtaining  this  result, 

we  assume  that  the  radiation  is  produced  via  a  Landau  resonance.  We 

also  assume,  for  this  calculation,  that  the  pitch  angle  of  the  electrons 

is  10°.  The  actual  pitch  angles  varied  from  0°  to  20°;  however,  the 

results  are  relatively  insensitive  to  pitch  angles  in  this  range.  From 

Figure  9  it  can  be  seen  that  the  most  intense  radiation  occurs  between 

the  electron  cyclotron  frequency  and  the  lower  hybrid  frequency,  f^HR* 

Outside  this  range  the  power  drops  by  a  factor  of  104 .  Note  that  this 

frequency  range  corresponds  rather  well  to  the  frequency  range  of  the 

radiation  observed  by  the  PDF.  Multiplying  the  power  from  each  electron 

by  the  number  of  electrons  in  the  first  200  meter  segment  of  the  beam 

(3  x  10^2  particles)  yields  (4-r’)  ~  lO-'^  W/Hz  in  the  frequency  range 

al  total 


from  fc  to  f^HR*  These  power  spectral  densities  are  much  lower  than 
the  measured  power  spectral  densities,  by  about  a  factor  of  107 
(compare  with  Figure  6,  where  dP/df  ~  10-9  W/Hz).  Therefore,  an  inco¬ 
herent  process  cannot  account  for  the  measured  wave  powers.  Some 
coherent  wave  process  must  be  involved.  In  Taylor  and  Shawhan's 
[1973]  analyses  of  the  generation  of  VLF  hiss  emissions  by  auroral 
electron  beams,  the  calculated  powers  for  the  incoherent  Cerenkov 
process  were  found  to  be  a  factor  of  102  -  103  lower  than  those 


mvm 
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COHERENT  GENERATION  MECHANISMS 


As  concluded  in  the  previous  section,  some  coherent  process  must 
be  involved  in  the  whistler-mode  wave  generation  from  the  SL-2  electron 
beam.  Coherent  processes  can  be  divided  into  two  classes:  direct  and 
indirect.  Direct  mechanisms  involve  the  direct  conversion  of  energy 
from  an  unstable  particle  distribution  to  electromagnetic  radiation; 
whereas  indirect  mechanisms  involve  the  intermediate  generation  of  one 
or  more  electrostatic  modes  which  are  coupled  to  the  escaping  electro¬ 
magnetic  radiation.  This  section  will  discuss  possible  direct  and 
indirect  mechanisms  that  may  explain  the  whistler-mode  radiation. 

Since  an  unstable  electron  distribution  is  present  in  the  beam  the 
escaping  electromagnetic  radiation  may  result  from  direct  conversion  of 
the  beam  energy  to  electromagnetic  radiation.  Such  a  mechanism  has 
been  proposed  by  Maggs  [1976]  for  the  generation  of  auroral  hiss.  In 
his  model,  incoherent  Cerenkov  radiation  produced  by  an  auroral  elec¬ 
tron  beam  is  directly  amplified  via  a  whistler-mode  plasma  instability 
within  the  beam.  It  seems  reasonable  that  this  wave  generation  mecha¬ 
nism  could  be  applied  to  the  whistler-mode  waves  emitted  from  the  SL-2 
electron  beam;  however,  a  problem  arises  in  doing  so.  Unlike  auroral 
beams,  the  path  length  for  wave  growth  in  the  SL-2  beam  is  very  short, 
only  two  to  three  electron  cyclotron  radii  (6  to  9  meters).  Using  the 
Landau  resonance  condition  and  the  fact  that  the  emission  is 
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propagating  near  the  resonance  cone,  we  can  derive  an  expression  for 
the  wavelength  of  the  whistler-mode  radiation, 


X 


T"  cos 
f  Res 


(12) 


which,  for  the  nominal  parameters  has  a  value  of  about  20  meters.  This 
wavelength  is  greater  than  the  path  length,  which  completely  invalidates 
any  mechanism  involving  exponential  growth.  Even  if  that  were  not  the 
case,  for  typical  whistler-mode  group  velocities  of  107  m/sec,  the  amount 
of  time  the  wave  spends  in  the  beam  is  so  short,  only  about  10-6  sec, 
that  unreasonably  high  growth  rates  (y  >  u>c  =  106  sec-1)  would  be 
required  to  generate  the  radiation.  No  whistler-mode  instability  is 
known  that  can  produce  such  large  growth  rates  from  realistic  electron 
distribution  functions.  These  same  conclusions  were  also  reached  by 
Jones  and  Kellogg  [1973]  in  their  paper  addressing  the  growth  rates  of 
whistler-mode  radiation  from  artificially-created  electron  beams. 

Next  we  consider  mechanisms  involving  the  intermediate  generation  of 
electrostatic  waves  in  the  beam.  Any  density  perturbation  or  bunch 
created  by  an  electrostatic  wave  in  the  beam  is  capable  of  emitting 
coherent  Cerenkov  radiation.  The  radiated  power  from  this  process  will 
have  a  frequency  spectrum  similar  to  that  of  a  single  radiating  electron; 
however,  the  wave  power  will  be  much  greater  since  the  emitted  power  goes 
as  N2  ,  where  N  is  the  number  of  electrons  in  a  bunch.  Coherent  Cerenkov 
radiation  from  a  bunched  beam  has  been  considered  previously  by  Bell 
[ 1968] . 


Beam-plasma  instabilities  are  known  to  be  capable  of  creating 

intense  electrostatic  waves  and  density  perturbations  in  the  beam.  We 

estimate  the  number  of  coherently  bunched  electrons  required  to  account 

for  the  observed  whistler-mode  radiation.  The  total  power  emitted  from 

the  electron  bunches  in  the  beam  is  (4r)  =  (— )  (AN)2a,  where 

df  total  df  e 
dP 

is  the  power  radiated  by  each  electron,  AN  is  the  typical  number 

of  electrons  in  a  bunch,  and  a  is  the  number  of  bunches  in  the  200 

meter  segment  of  the  beam.  Beam-plasma  instabilities  are  known  to 

create  an  electrostatic  wave  near  the  local  electron  plasma  frequency 

(3  MHz).  Such  an  emission  is,  in  fact,  observed  near  3  MHz  [see 

Gurnett  et  al.,  1986].  The  corresponding  wavelength  is  V^/fp  -  7 

meters,  which  is  assumed  to  be  the  approximate  length  of  each  bunch. 

This  wavelength  can  then  be  used  to  calculate  a,  the  number  of  bunches 

in  the  first  200  meters  of  the  beam.  This  number  is  a  =  29.  Since  the 

dP 

radiated  power  from  the  200  meter  beam  segment,  (t?)  ,  is  about 

df  total 

10~9  W/Hz,  is  about  10-29  W/Hz  and  a  -  29 ,  each  bunch  must 

ar  e 

contain  about  AN  =  2  x  109  electrons  in  order  to  account  for  the 
observed  radiated  power. 

Next  we  estimate  the  required  electric  field  strength  for  the 
electrostatic  wave  in  the  beam  that  forms  the  bunches.  Assuming  that 
the  beam  diameter  is  about  2  cyclotron  radii,  the  electron  number 
density  in  the  bunch  can  be  estimated  using  the  formula: 
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where  AL  is  the  bunch  length  and  rc  is  the  cyclotron  radius  (2  to  3 
meters).  The  required  number  density  is  found  to  be  about  An  =  1  x 
107  electrons/m3.  Note  that  the  fractional  density  perturbation  in 
the  beam  An/n0  is  only  about  0.01;  thus,  a  relatively  small  perturba¬ 
tion  can  account  for  the  measured  whistler-mode  power.  Poisson's 
equation  can  be  used  to  determine  the  magnitude  of  the  self-consistent 
electric  field  needed  to  generate  this  density  perturbation 

AE  _  eAn  ... . 

AL  *  (14) 

o 

From  (14),  we  estimate  that  an  electric  field  on  the  order  of  1-2  V/m 
is  needed  to  create  the  required  coherence  in  the  beam  electrons. 

Although  the  PDP  did  not  fly  directly  through  the  beam  during 
free  flight,  when  the  PDP  was  on  the  remote  manipulator  arm,  it  did 
provide  measurements  in  the  beam.  During  these  times,  intense 
field-aligned  electric-field  signals  near  fpe  were  measured  with 
amplitudes  greater  than  0.3  V/m,  sufficiently  large  to  saturate  the 
receiver.  This  value  is  within  a  factor  of  10  of  the  required 
amplitudes  needed  for  radiative  coherence  of  the  beam  electons.  The 
good  agreement  between  the  calculated  and  measured  electrostatic  field 
strengths  strongly  suggests  that  electron  hunches  generated  by  a 
beam-plasma  instability  can  account  for  the  observed  whistler-mode 


power. 
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CONCLUSIONS 


Like  auroral  hiss,  the  whistler-mode  radiation  from  the  SL-2  elec¬ 
tron  beam  was  found  to  be  progagating  near  the  resonance  cone  in  the 
same  direction  as  the  beam.  This  conclusion  was  reached  by  comparing 
the  measured  and  modeled  electric-f ield  polarizations  in  the  PDP  spin 
plane.  The  total  whistler-mode  power  radiated  by  the  beam  between  the 
shuttle  and  the  PDP  was  estimated  to  be  1.6  mW.  This  represents  a  con¬ 
version  of  about  3.2  x  10~5  of  the  beam  power  to  whistler-mode  radia¬ 
tion  in  a  path  length  of  200  meters.  The  linear  emissivity  was  also 
observed  to  decreased  with  increasing  distance  along  the  beam,  thereby 
indicating  that  most  of  the  whistler-mode  radiation  is  generated  rela¬ 
tively  close  to  the  source  of  the  beam. 

The  mechanism  that  best  explains  the  radiated  power  and  frequency 
range  of  the  whistler-mode  radiation  is  Cerenkov  radiation  from  coher¬ 
ently  radiating  electron  bunches  in  the  beam.  The  bunches  are  most 
likely  created  by  an  electrostatic  instability  in  the  beam.  The  esti¬ 
mated  electric  field  amplitudes  required  to  create  the  bunching  are 
consistent  with  the  field  strength  of  electrostatic  noise  measured  in 
the  beam. 

Future  work  includes  calculating  the  radiated  power  from  a  model 
electron  beam  obtained  from  a  particle  simulation.  The  calculated 
power  can  then  be  compared  to  the  measured  whistler-mode  power  to  con- 
f i rm  that  coherent  radiation  from  electron  bunches  generate  the  signal 
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FIGURE  CAPTIONS 


Figure  1  A  frequency  vs.  time  spectogram  from  the  POP  plasma 

wave  instrument  showing  intense  emissions  during  a  D.C. 
electron  gun  firing.  The  funnel-shaped  structure  that 
extends  from  the  electron  cyclotron  frequency,  fc,  to 
about  30  MHz  is  whistler-mode  radiation  from  the  beam. 

Figure  2  This  diagram  shows  the  index  of  refraction  surface  for 

the  whistler  mode  and  the  associated  E,  k,  and  Vg  vectors 
for  propagation  near  the  resonance  cone  (f)  -  ARes).  For 
propagation  near  the  resonance  cone,  k  and  E  are  parallel 
and  nearly  perpendicular  to  vg.  In  this  limit  E  is 
linearly  polarized  and  quasi-electrostatic. 

Figure  3  This  figure  shows  the  whistler-mode  ray  paths  from  a 

point  source  radiator.  Note  that  a  passing  spacecraft 
will  detect  radiation  near  the  gyrof requency ,  first,  then 
detect  lower  frequency  emissions  as  it  nears  the  source. 

A  well-defined  frequency  versus  time  funnel-shaped 
spectra  results  from  this  spacecraf t/source  encounter. 

Figure  4  This  diagram  shows  the  ray  path  and  E,  k,  and  Vg  vectors 
used  to  confirm  the  electric  field  polarization.  The 
assumed  electric  field  is  projected  into  the  PDP  spin 
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plane  and  the  angle  relative  to  the  projection  of  the  sun 
vector  is  calculated.  The  projected  electric  field 
direction  can  then  be  compared  to  the  measured  directions 
calculated  from  spin  modulation  maxiraums  in  the  electric 
field  intensity  (see  Figure  4). 

Figure  5(a),  (b),  (c),  and  (d)  These  plots  show  the  relative  directions 
of  the  computed  and  measured  electric-field  vectors  in  the 
PDP  spin  plane  for  the  562  kHz,  311  kHz,  178  kHz,  and  100 
kHz  frequency  channels.  The  dots  represent  the  computed 
electric  field  directions  assuming  that  the  wave  vector  is 
near  the  resonance  cone  with  k*v{,  >  0,  and  the  x's  repre¬ 
sent  measured  electric-field  directions.  The  close  agree¬ 
ment  between  the  measured  and  modeled  directions  indicates 
that  the  whistler-mode  radiation  is  propagating  near  the 
resonance  cone  in  the  same  direction  as  the  beam. 

Figure  6  This  diagram  shows  the  integration  surface  used  to 

calculate  the  power  emitted  from  the  beam  in  the  whistler 
mode.  At  closest  approach,  the  PDP  passed  within  3  meters 
of  the  beam  at  a  distance  of  about  200  meters  from  the 
shuttle. 

Figure  7  The  calculated  power  spectral  density  from  the  beam  in  the 
whistler  mode  is  shown  as  a  function  of  frequency. 

Figure  8(a)  and  (b)  The  linear  emissivity,  dP/dfdi,  is  shown  as  a 
function  of  the  distance,  L,  along  the  beam  for  the  562 
kHz  and  311  kHz  frequency  channels.  Note  that  the 
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emissivity  starts  to  decrease  rapidly  beyond  about  IDO 
meters. 

Figure  9  The  power  spectra  from  a  single  electron  radiating  via 
the  Cerenkov  processes  is  shown  in  a  plasma  environment 
similar  to  that  surrounding  the  SL-2  beam.  These 
calculations  assume  the  wave/beam  interaction  is  by  a 
Landau  resonance  process  and  that  the  particle  pitch 
angle  is  10°.  This  power  calculation  is  based  on 
formulas  derived  by  Mansfield  [1967]. 
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